------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
 opened on:  23 Nov 2021, 16:54:58

. log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
 paused on:  23 Nov 2021, 16:54:58
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
resumed on:  23 Nov 2021, 16:54:58

.   oprobit cat bi fc sc tn tp turb id1 id2 id3 id4 id5 [fw=ec], vce(cluster id) ;

note: id5 omitted because of collinearity
Iteration 0:   log pseudolikelihood = -52422.971  
Iteration 1:   log pseudolikelihood = -40187.476  
Iteration 2:   log pseudolikelihood = -40062.411  
Iteration 3:   log pseudolikelihood = -40062.104  
Iteration 4:   log pseudolikelihood = -40062.104  

Ordered probit regression                       Number of obs     =     50,000
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -40062.104               Pseudo R2         =     0.2358

                                     (Std. Err. adjusted for 5 clusters in id)
------------------------------------------------------------------------------
             |               Robust
         cat |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          bi |   .9109499   .1581352     5.76   0.000     .6010106    1.220889
          fc |   .0000686   .0000262     2.62   0.009     .0000172    .0001199
          sc |   .0011447   .0001006    11.37   0.000     .0009475     .001342
          tn |   .1027869   .0170235     6.04   0.000     .0694215    .1361523
          tp |   .2245749   .0987958     2.27   0.023     .0309386    .4182112
        turb |   .0081166   .0032368     2.51   0.012     .0017726    .0144606
         id1 |  -.0041305   .0023902    -1.73   0.084    -.0088152    .0005542
         id2 |    .177707   .0066523    26.71   0.000     .1646688    .1907453
         id3 |  -.6689765   .0154266   -43.37   0.000    -.6992119    -.638741
         id4 |  -.6708903   .0153217   -43.79   0.000    -.7009204   -.6408603
         id5 |          0  (omitted)
-------------+----------------------------------------------------------------
       /cut1 |    4.86002    .815709                       3.26126    6.458781
       /cut2 |   6.587022   .8098423                      4.999761    8.174284
------------------------------------------------------------------------------

.   log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
 paused on:  23 Nov 2021, 16:54:59
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
resumed on:  23 Nov 2021, 16:54:59

.   oprobit cat bi fc sc tn tp turb id1 id2 id3 id4 id5 [fw=mw], vce(cluster id) ;

note: id5 omitted because of collinearity
Iteration 0:   log pseudolikelihood = -52084.742  
Iteration 1:   log pseudolikelihood = -41695.858  
Iteration 2:   log pseudolikelihood = -41636.127  
Iteration 3:   log pseudolikelihood = -41635.979  
Iteration 4:   log pseudolikelihood = -41635.979  

Ordered probit regression                       Number of obs     =     50,000
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -41635.979               Pseudo R2         =     0.2006

                                     (Std. Err. adjusted for 5 clusters in id)
------------------------------------------------------------------------------
             |               Robust
         cat |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          bi |   .3167587   .1037627     3.05   0.002     .1133876    .5201299
          fc |   .0000479   .0000167     2.87   0.004     .0000152    .0000807
          sc |   .0002734   .0005562     0.49   0.623    -.0008168    .0013636
          tn |   .0741036   .0242275     3.06   0.002     .0266185    .1215886
          tp |   .4031254    .099961     4.03   0.000     .2072055    .5990453
        turb |   .0232159   .0043526     5.33   0.000      .014685    .0317467
         id1 |  -.9240597   .0290829   -31.77   0.000    -.9810611   -.8670582
         id2 |  -.6690151   .0216566   -30.89   0.000    -.7114612    -.626569
         id3 |  -.9172759   .0289536   -31.68   0.000    -.9740238   -.8605279
         id4 |  -.8632061   .0263951   -32.70   0.000    -.9149395   -.8114728
         id5 |          0  (omitted)
-------------+----------------------------------------------------------------
       /cut1 |   1.895651   .6842908                      .5544659    3.236837
       /cut2 |   3.422477   .6126184                      2.221767    4.623187
------------------------------------------------------------------------------

.   log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
 paused on:  23 Nov 2021, 16:54:59
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
resumed on:  23 Nov 2021, 16:54:59

.   oprobit cat bi  fc  sc  tn  tp  turb  id6 id7 id8 [fw=hr], vce(cluster id) ;

note: id8 omitted because of collinearity
Iteration 0:   log pseudolikelihood = -32642.626  
Iteration 1:   log pseudolikelihood = -25193.765  
Iteration 2:   log pseudolikelihood = -25131.515  
Iteration 3:   log pseudolikelihood = -25131.389  
Iteration 4:   log pseudolikelihood = -25131.389  

Ordered probit regression                       Number of obs     =     30,000
                                                Wald chi2(0)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -25131.389               Pseudo R2         =     0.2301

                                     (Std. Err. adjusted for 3 clusters in id)
------------------------------------------------------------------------------
             |               Robust
         cat |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          bi |   .2375807   .0599486     3.96   0.000     .1200835    .3550779
          fc |   .0002518   .0000937     2.69   0.007     .0000682    .0004354
          sc |   .0021633   .0003986     5.43   0.000     .0013821    .0029446
          tn |   .0695236   .0132899     5.23   0.000     .0434759    .0955713
          tp |  -.0760291   .1052173    -0.72   0.470    -.2822513    .1301931
        turb |   .0130247   .0055485     2.35   0.019     .0021498    .0238996
         id6 |  -.9751387   .0627128   -15.55   0.000    -1.098054   -.8522238
         id7 |   .1404157   .0355237     3.95   0.000     .0707906    .2100408
         id8 |          0  (omitted)
-------------+----------------------------------------------------------------
       /cut1 |   1.970425   .6804315                      .6368041    3.304047
       /cut2 |   3.258697    .654771                       1.97537    4.542025
------------------------------------------------------------------------------

.     log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
 paused on:  23 Nov 2021, 16:54:59
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
resumed on:  23 Nov 2021, 16:54:59

. logit ce cost          ec_g ec_p hr_g hr_p md_g md_p age gender white college [pw=wgt_final], cluster(resp_id) ;

Iteration 0:   log pseudolikelihood = -4191743.7  
Iteration 1:   log pseudolikelihood = -3687122.1  
Iteration 2:   log pseudolikelihood = -3686029.4  
Iteration 3:   log pseudolikelihood = -3686029.3  

Logistic regression                             Number of obs     =      7,584
                                                Wald chi2(11)     =     512.68
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -3686029.3               Pseudo R2         =     0.1206

                            (Std. Err. adjusted for 1,896 clusters in resp_id)
------------------------------------------------------------------------------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        cost |  -.0561872   .0034198   -16.43   0.000    -.0628898   -.0494846
        ec_g |   .0162029    .003597     4.50   0.000      .009153    .0232528
        ec_p |  -.0152104   .0038619    -3.94   0.000    -.0227796   -.0076411
        hr_g |   .0113868   .0029719     3.83   0.000      .005562    .0172116
        hr_p |  -.0248146   .0087398    -2.84   0.005    -.0419444   -.0076848
        md_g |   .0078538   .0022209     3.54   0.000      .003501    .0122066
        md_p |  -.0021807   .0047154    -0.46   0.644    -.0114227    .0070614
         age |  -.0162309    .002857    -5.68   0.000    -.0218305   -.0106312
      gender |  -.1008659   .1026613    -0.98   0.326    -.3020784    .1003466
       white |  -.2496793   .1215481    -2.05   0.040    -.4879092   -.0114494
     college |   .2653966   .1025011     2.59   0.010     .0644982     .466295
       _cons |   1.288001   .2250702     5.72   0.000     .8468716    1.729131
------------------------------------------------------------------------------

. nlcom wtp : -12*(_b[ec_g]-_b[ec_p]+_b[hr_g]-_b[hr_p]+_b[md_g]-_b[md_p])/_b[cost] ;

         wtp:  -12*(_b[ec_g]-_b[ec_p]+_b[hr_g]-_b[hr_p]+_b[md_g]-_b[md_p])/_b[cost]

------------------------------------------------------------------------------
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         wtp |   16.58366   2.791154     5.94   0.000      11.1131    22.05422
------------------------------------------------------------------------------

. logit ce cost cost_dis ec_g ec_p hr_g hr_p md_g md_p age gender white college [pw=wgt_final], cluster(resp_id) ;

Iteration 0:   log pseudolikelihood = -4191743.7  
Iteration 1:   log pseudolikelihood = -3684516.7  
Iteration 2:   log pseudolikelihood =   -3683448  
Iteration 3:   log pseudolikelihood = -3683447.9  

Logistic regression                             Number of obs     =      7,584
                                                Wald chi2(12)     =     514.22
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -3683447.9               Pseudo R2         =     0.1213

                            (Std. Err. adjusted for 1,896 clusters in resp_id)
------------------------------------------------------------------------------
             |               Robust
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        cost |  -.0610311   .0046437   -13.14   0.000    -.0701326   -.0519295
    cost_dis |   .0401019   .0250425     1.60   0.109    -.0089804    .0891843
        ec_g |   .0162943   .0035968     4.53   0.000     .0092448    .0233438
        ec_p |  -.0151874   .0038603    -3.93   0.000    -.0227535   -.0076213
        hr_g |   .0114533   .0029718     3.85   0.000     .0056287    .0172779
        hr_p |  -.0247109   .0087397    -2.83   0.005    -.0418405   -.0075814
        md_g |   .0079231   .0022232     3.56   0.000     .0035657    .0122805
        md_p |  -.0021735   .0047079    -0.46   0.644    -.0114009    .0070539
         age |  -.0162039   .0028559    -5.67   0.000    -.0218013   -.0106064
      gender |  -.0959678   .1028224    -0.93   0.351     -.297496    .1055604
       white |   -.257121   .1215866    -2.11   0.034    -.4954263   -.0188157
     college |   .2558948   .1027845     2.49   0.013     .0544409    .4573488
       _cons |   1.293062   .2249336     5.75   0.000     .8522005    1.733924
------------------------------------------------------------------------------

. nlcom wtp : -12*(_b[ec_g]-_b[ec_p]+_b[hr_g]-_b[hr_p]+_b[md_g]-_b[md_p])/(_b[cost]+_b[cost_dis]/(1+ 9.584872)) ;

         wtp:  -12*(_b[ec_g]-_b[ec_p]+_b[hr_g]-_b[hr_p]+_b[md_g]-_b[md_p])/(_b[cost]+_b[cost_dis]/(1+ 9.584872))

------------------------------------------------------------------------------
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         wtp |   16.29752   2.734841     5.96   0.000     10.93733    21.65771
------------------------------------------------------------------------------

.   /* WTP FOR PERSON LIVING AVERAGE DISTANCE FROM CENTROID CLEANED-UP STREAMS */
> nlcom wtp : -12*(_b[ec_g]-_b[ec_p]+_b[hr_g]-_b[hr_p]+_b[md_g]-_b[md_p])/(_b[cost]+_b[cost_dis]/(1+ 1)) ;

         wtp:  -12*(_b[ec_g]-_b[ec_p]+_b[hr_g]-_b[hr_p]+_b[md_g]-_b[md_p])/(_b[cost]+_b[cost_dis]/(1+ 1))

------------------------------------------------------------------------------
          ce |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         wtp |   22.76496   6.681436     3.41   0.001     9.669586    35.86034
------------------------------------------------------------------------------

.          /* WTP FOR PERSON LIVING ONE MILE         FROM CENTROID CLEANED-UP STREAMS */
> log off ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
 paused on:  23 Nov 2021, 16:55:00
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
resumed on:  23 Nov 2021, 16:55:00

. sum wtp [w=wgt_final] if            dis < 1  & stream == "Crabtree & Walnut" & scenario == 1 ;
(analytic weights assumed)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
         wtp |       9  4349.23466    173.3154    21.9099   154.8881   216.1068

. sum wtp [w=wgt_final] if dis >= 1 & dis < 2  & stream == "Crabtree & Walnut" & scenario == 1 ;
(analytic weights assumed)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
         wtp |      17  10358.9902    137.2341   6.194329    130.292   148.5382

. sum wtp [w=wgt_final] if dis >= 2 & dis < 5  & stream == "Crabtree & Walnut" & scenario == 1 ;
(analytic weights assumed)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
         wtp |     157  117815.364    117.6426   3.789507    113.424   129.0086

. sum wtp [w=wgt_final] if dis >= 5 & dis < 10 & stream == "Crabtree & Walnut" & scenario == 1 ;
(analytic weights assumed)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
         wtp |     306  246365.809    109.5802    1.68179   107.4072   113.3926

. sum wtp [w=wgt_final] if dis >=10 & dis < 30 & stream == "Crabtree & Walnut" & scenario == 1 ;
(analytic weights assumed)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
         wtp |     296  263149.673    105.8485   .9240938   103.8387   107.3834

. sum dis, detail ;

                 Distance to Impacted Creeks
-------------------------------------------------------------
      Percentiles      Smallest
 1%     .7706293       .2335068
 5%     2.452453       .2335068
10%     3.380927       .2335068       Obs               3,140
25%     5.264293       .2335068       Sum of Wgt.       3,140

50%     8.419018                      Mean           8.933176
                        Largest       Std. Dev.      4.567334
75%     11.88053       22.94525
90%     15.30865       22.94525       Variance       20.86054
95%     17.10397       22.94525       Skewness       .4591577
99%     21.02459       22.94525       Kurtosis       2.661109

. sort stream scenario ;

. di "HOUSEHOLD-LEVEL WAKE COUNTY BENEFITS" ;
HOUSEHOLD-LEVEL WAKE COUNTY BENEFITS

. collapse (mean) wtp [w=wgt_final], by(stream scenario) ;
(analytic weights assumed)

. list ;

     +------------------------------------------+
     |            stream   scenario         wtp |
     |------------------------------------------|
  1. | Crabtree & Walnut          1   110.40807 |
  2. | Crabtree & Walnut          2   33.269011 |
  3. |    Middle & Swift          1   75.444448 |
  4. |    Middle & Swift          2   26.933905 |
     +------------------------------------------+

. replace wtp = wtp*400172 ;
(4 real changes made)

.    /* WAKE COUNTY HOUSEHOLDS IS 400,172 - https://www.census.gov/quickfacts/fact/table/wakecountynorthcarolina/PST045219 */
> di "AGGREGATE WAKE COUNTY BENEFITS" ;
AGGREGATE WAKE COUNTY BENEFITS

. list ;

     +-----------------------------------------+
     |            stream   scenario        wtp |
     |-----------------------------------------|
  1. | Crabtree & Walnut          1   44182217 |
  2. | Crabtree & Walnut          2   13313327 |
  3. |    Middle & Swift          1   30190756 |
  4. |    Middle & Swift          2   10778195 |
     +-----------------------------------------+

. log close ;
      name:  <unnamed>
       log:  C:\Users\rhhaefen\Dropbox\EPA_WQ_Grant_Survey_Files\Final_Data_Public\log\4-case_study.log
  log type:  text
 closed on:  23 Nov 2021, 16:55:00
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
